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ABSTRACT 


Wind  tunnel  studies  of  particulate  lofting  were  conducted 
at  free  stream  speeds  of  between  112  and  376  fps.  Measurements 
were  made  of  the  rate  at  which  particulates  were  lofted ,  of  the  wind 
velocity  profile  and  of  the  number  density  of  airborne  particulates. 
The  measurement  bearing  the  most  practical  significance  was  that 
of  the  lofting  rate;  the  measured  lofting  rate  was  only  1/10  that 
of  the  maximum  rate  predicted  by  current  theory,  although  the 
lofting  rate  was  found  to  be  directly  proportional  to  the  shear 
stress  velocity  in  accord  with  the  theory.  The  velocity  profile 
within  the  boundary  layer  but  at  heights  greater  than  the  effective 
roughness  created  by  the  particulate  flow  is  of  the  form  given  by 
Prandtl's  law  for  flow  past  a  wall.  Boundary  layer  calculations 
in  conjunction  with  the  measurements  show  that  Owen's  form  for  the 
flow  law  cannot  hold  near  the  leading  edge  of  the  boundary  layer. 
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Defense  Atomic  Support  Agency  (DASA)  under  Contract  DASA01-70-C-0041 
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DASA  contracting  officer's  representative  was  Mr.  Clifford  B. 
McFarland. 

The  Principal  investigator  for  the  program  was 
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W.  Gordy  who  built  the  apparatus,  L.  Walitt,  M.  Chapman  and 
L.  Glenn  who  aided  in  the  aerodynamic  design  of  the  experiments 
and  the  staff  of  the  Fluid  Dynamics  Laboratory  of  the  Lockheed- 
California  Company  Rye  Canyon  Research  Laboratory  where  the 
experiments  were  performed. 
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NOMENCLATURE 


area 

area  at  sonic  flow 

local  friction  coefficient 

diameter 

acceleration  of  gravity 

mass  of  particulates  per  unit  area 

mass  lofting  rate  lbm/Cft^-sec) 

reciprocal  of  exponent  in  power  law  velocity 
profile 

vertical  velocity  of  particulate 

free  stream  velocity 

shear  stress  velocity 

axial  coordinate 

vertical  coordinate 

boundary  layer  thickness 

displacement  thickness 

momentum  thickness 

gas  density 

time  constant 

diameter  of  particulate 

light  intensity  -  watts/ steradian 

initial  light  intensity  -  watts/steradian 

total  scattering  coefficient  (total  scattering 
cross  section  -f  geometric  cross  section) 

path  length 

number  density  -  nuraber/cm3 
radius  of  particulate 
X/2nr 

wavelength 

angle 

angular  radius  of  Airy's  disc 
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SECTION  1 


INTRODUCTION  AND  SUMMARY 

In  the  past,  attempts  have  been  made  to  calculate  the 
mass  and  velocity  of  a)  debris  thrown  out  directly  from  the 
crater  resulting  from  a  nuclear  explosion  and  b)  the  debris 
lofted  outside  the  crater  by  the  high  speed  wind  produced  by 
the  same  nuclear  explosion.  The  latter  problem  is  the  concern 
of  this  research.  The  lofting  calculations  were  based  in  part 
upon  meteorological  and  laboratory  data  obtained  with  low  speed 
winds  which  were  extrapolated  to  the  very  high  wind  speeds 
characteristic  of  nuclear  explosion  phenomena,  and  in  part  upon 
a  theoretical  model  of  the  lofting  process.  The  purpose  of 
this  program  was  to  obtain  for  the  steady  state  lofting  process 
high  speed  laboratory  data  which  could  be  used  to  Increase  the 
accuracy  of  the  theoretical  model.  The  nuclear  problem  is 
very  complicated;  the  problem  addressed  here  was  highly  ideal¬ 
ized  In  an  attempt  to  uncover  and  understand  the  basic  aero¬ 
dynamics  of  the  lofting  process.  To  this  end  the  experiments 
were  performed  In  a  wind  tunnel  with  the  ground  surface 
idealized  by  a  bed  of  finely  graded,  noncohesive  silica  particles 
For  the  tests  to  have  meaning  In  the  context  of  lofting  on  a 
meteorological  scale  the  wind  tunnel  was  made  large  enough  so 
that  the  boundary  layer  did  not  occupy  an  appreciable  fraction 
of  the  tunnel  height  and  so  that  the  number  of  particulates 
striking  the  roof  of  the  tunnel  was  so  small  that  their 
presence  did  not  perturb  the  flow. 

Three  tests  were  performed  at  free  stream  speeds  of  111.8 
237.8  and  371.3  fps  which  placed  the  particulate  motion  in  the 
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suspension  regime  (the  ratio  of  the  drag  force  based  on  the  shear 
stress  velocity  to  the  weight  of  the  particulate  was  3,3,  17.4, 
and  88.0  at  the  three  wind  speeds).  Measurements  were  made  of  the 
rate  at  which  particulates  were  lofted,  of  the  velocity  profile, 
and  of  the  number  density  of  lofted  particulates.  The  measurement 
bearing  the  most  practical  significance  was  that  of  the  lofting 
rate;  the  measured  lofting  rate  was  only  1/10  the  maximum  rate  pre¬ 
dicted  by  the  theoretical  model.  The  lofting  rate  was  found  to  be 
directly  proportional  to  the  shear  stress  velocity,  in  accord  with 
the  model.  The  difference  between  the  calculated  and  the  measured 
lofting  rate  does  not  imply  that  the  model  is  invalid,  but  rather  that 
the  experiments  can  be  used  within  the  framework  of  the  model  to 
set  the  efficiency  of  conversion  of  the  horizontal  impulse  deliv¬ 
ered  to  the  ground  surface  by  the  wind  to  the  vertical  momentum  of 
particulates.  The  calculation,  aimed  at  maximizing  the  lofting 
rate,  set  this  conversion  efficiency  at  100%;  the  experiments  imply 
that  the  conversion  process  is  only  10%  efficient.  Whether  or 
not  all  of  the  details  of  the  model  are  correct  requires  more 
extensive  experimentation.  The  design  phase  of  the  present  work, 
however  did  point  out  a  deficiency  in  the  calculation  presented  in 
Reference  8;  namely,  Owen's  relation,  Eq.  4,  was  applied  without 
regard  to  the  actual  height  of  the  boundary  layer.  The  calcula¬ 
tion  of  gas  motion,  which  was  inviscid,  did  not  account  for  the 
growth  of  the  boundary  layer  behind  the  initial  blast  wave  and  the 
subsequent  velocity  profile  arising  from  the  afterwinds. 

The  measured  velocity  profiles  were  interpreted  in  accord 
with  the  law  for  flow  past  a  wall,  with  Owen's*  relation  providing 
a  bridge  between  the  test  at  111.8  fps  and  the  low  speed  tests  per¬ 
formed  in  the  saltation  regime  by  other  investigators.  Owen's  rela¬ 
tion  did  not  hold  for  the  two  higher  speed  tests;  boundary  layer 
theory  indicated  that  Owen's  law  cannot  apply  near  the  leading  edge 
of  the  boundary  layer. 


2 


SECTION  2 


EXPERIMENTAL  TECHNIQUES 


2.1  TEST  SETUP 

Steady  state  lofting  tests  were  performed  inside  a  rectan¬ 
gular  cross  section  wind  tunnel  with  a  40"  height,  IP  '  width  and 
26' 8"  length  (Figure  1).  A  4' “long  circular- to- rectangular  transi¬ 
tion  section  joined  the  tunnel  to  Lockheed-California  Company's 
48"  cross  wind  nozzle.  A  particulate  bed  approximately  4"  deep  at 
the  start  of  each  test  was  composed  of  AFS  50-70  Ottawa  silica 
testing  sand  (mean  particulate  diameter:  .25  mm).  Design  details 
of  the  steady  state  tunnel  along  with  that  of  a  transient  tunnel 
are  given  in  Appendix  I.  Lofting  measurements  were  made  within  a 
diagnostics  section  which  extended  between  16'  and  18' 6"  from  the 
leading  edge  of  the  particulate  bed.  The  location  of  the  diagnostics 
section  was  determined  as  a  compromise  between  placing  the  measure¬ 
ments  station  at  least  10  hydraulic  diameters  downstream  of  the 
entrance  and  fitting  the  tunnel  into  the  available  space.  The  veloc¬ 
ity  profile  of  the  wind  was  obtained  from  differential  pressure 
measurements  utilizing  a  rake  containing  8  coaxial  pitot  tubes  and  3 
static  tubes.  The  density  profile  of  the  lofted  particulates  was 
inferred  from  light  extinction  measurements.  The  position  of  the 
eroding  ground  surface  needed  for  the  velocity  and  density  measure¬ 
ments  was  obtained  photographically  using  a  Locara  camera.  The  loft¬ 
ing  rate  was  obtained  from  measurements  of  the  change  in  the  level 
of  the  particulate  bed  during  a  test  and  the  test  time  (Appendix  II). 
Pressure,  light  intensity  and  gas  temperature  measurements  were 
recorded  using  Lockheed-California  Company's  central  data  system. 

Each  channel  employed  a  10  cps  low  pass  filter  and  was  sampled  at 
a  rate  of  40/sec. 
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2.2 


DIAGNOSTICS 


2.2.1  Particulate  Number  Density 

Under  certain  conditions,  which  are  described  below, 
the  density  of  particulates  contained  within  a  beam  of  radia¬ 
tion  can  be  ascertained  by  application  of  the  Bouguer-Beer 
relation, 


I  *  I  exp(-niTr2K  l)  (1) 

O  X 

to  the  results  of  light  transmission  experiments.  To  do  so 
requires,  first  of  all,  that  the  Bouguer-Beer  relation  hold 
for  the  experimental  setup.  The  main  effort  in  this  phase  of 
the  research  was  to  design  the  extinction  experiment  such  that 
Eq.  1  could  be  applied  in  a  meaningful  way  to  the  results  of 
the  experiment. 

It  is  clear  from  the  derivation  of  Eq.  1  that  in  each 
differential  element  of  length  dx  a  quantity  of  light  equal  to 
InTnraK  dx  is  scattered  out  of  the  beam,  but  that  no  radiation 
is  scattered  into  the  beam;  hence,  strict  application  of  this 
relation  requires  that  the  particulates  be  sufficiently  dif¬ 
fuse  or  that  the  path  length  be  sufficiently  short  that  multiple 
scattering  effects  be  negligible.  A  light  transmitter  and  re¬ 
ceiver  were  built  with  a  light  path  of  variable  length  in  order 
to  detect  multiple  scattering  effects.  It  is  assumed  that 
if  multiple  scattering  is  involved  in  the  scattering  process 
when  the  beam  length  is  equal,  say,  to  the  full  width  of  the 
tunnel,  then  there  exists  some  shorter  distance  for  which  only 
single  scattering  is  involved.  Secondly,  the  total  scattering 

coefficient,  K  ,  must  be  known.  In  general,  K  depends  on  both 

x 
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the  ratio  (denoted  by  o/)  of  the  particle  circumference  to  the 
wavelength  of  the  radiation,  and  on  the  refractive  index  of  the 
particle  relative  to  that  of  the  surrounding  medium.  The 
refractive  index  will  be  a  real,  quantity  for  a  lossless 
material,  but  will  be  complex  for  a  lossy  (absorbing)  material. 
For  a  lossless  material,  Kx  is  an  oscillatory  function  of  a 
which  approaches  the  value  2  for  a»  1,  while  for  a  lossy 
material,  K  is  heavily  damped  but  still  approaches  K  *  2  for 
large  a.  For  the  experiments  at  hand,  typically  r  =  125m-  and 
taking  X  =  550  nm  (green  light)  as  an  average  value  yields 
a  =  1430.  For  this  large  value  of  a,  K  -  2  independent  of  ct 

X 

and  the  refractive  index.  That  Kx  approaches  2  rather  than  1 
can  be  seen  more  readily  by  application  of  geometrical  and 
physical  optics  to  the  scattering  process,  than  by  applying  the 
rigorous  Mie  theory  for  a  »  1.  For  transparent  spherical 
particles  the  scattered  light  is  composed  of  reflected, 
refracted,  and  diffracted  light.  The  reflected  and  refracted 
components  result  from  light  incident  on  the  projected  area  of 
the  particle  and  would  yield  K  -  1  if  these  were  the  only 
components  of  the  scattered  light.  However,  by  Babinet's 
principle  an  amount  of  light  equal  to  that  incident  on  an 
opaque  disc  with  area  equal  to  that  of  the  projected  area  of 
on  opaque  particle  is  diffracted  around  the  particle.  Thus, 

K  -  2.  Exact  Mie  scattering  calculations  have  shown  that 

A 

diffraction  theory  for  an  opaque  circular  disc  is  a  very  good 

2 

approximation  for  diffraction  by  a  transparent  particle. 
Approximately  827.  of  the  diffracted  light  falls  within  Airy's 
disc.  A  good  rule  of  thumb,  then,  for  practical  design  is  that 
the  diffracted  light  falls  within  Airy's  disc.  A  more  exact 
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3 

relation  based  on  Rayleigh's  result  for  the  fraction  of  light 
falling  outside  a  cone  of  half  angle  cp  gives  an  apparent  total 
scattering  coefficient  of  1  +  C *Jo  (acp)  ]s  +  [Ji(a<p)]s  for  cp  small 
(sin  cp  ~  cp) ,  where  Jq  and  J,.  are  Bessel  functions  of  the  first 
kind,  of  order  zero  and  one,  respectively.  The  angular  radius 
of  Airy's  disc,  given  by  cp^  =  1.22  X/d,  is  0.152  degrees  for  a 
250m,  particle  illuminated  with  green  light  (X  »  550  run) .  Unless 
the  angular  aperture  of  the  light  detector  is  considerably 
smaller  than  cp^  the  apparent  total  scattering  coefficient  will 
be  less  than  2.  The  approach  adopted  here  is  to  make  the 
angular  aperture  of  the  receiver  optics  several  times  that  of 
cp^  to  ensure  that  all  of  the  diffracted  light  is  collected. 

The  angular  aperture  must  not  be  so  large,  however,  as  to 
collect  a  significant  amount  of  the  refracted  and  reflected 
light.  For  a  spherical  particle  geometric  optics  implies 
that  the  ratio  of  the  intensities  of  refracted  and  incident  light 

4 

varies  with  angle  according  to  the  formula 


.k  -  (1  -  ki )8  sin  2iT  di1 
IQ  2  sin  cp  dtp 


where  cp  »  2(i*  -  rx) 


ki 


1  [sinflk-L-Eil  +  -  ri>~ 

2  |_sin  (ix  +  r*)  tana(ii  +  rx). 
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rT  is  the  angle  of  refraction  and  ix  is  the  angle  of  incidence. 
The  angular  variation  of  the  reflected  light  is 


where 


■f*  -  ik, 

O 

*  I  rslna(ia  -  ra)  tan^^-^l 
2  [sina(ia  +  ra)  tan2  (is  +  rs)J 

=  rr  -  2i?  ; 


(3) 


r2  is  the  angle  of  reflection  corresponding  to  i8 . 

The  dependence  on  angle  of  the  ratio  of  the  sum  of  the  refracted 
and  reflected  intensities  to  the  incident  intensity  (Figure  2), 
indicates  that  the  measurement  error  will  be  small  provided  that 
the  angular  aperture  is  less  than  1  degree.  It  has  already  been 

shown  that  an  angle  of  0.152  degrees  is  sufficient  to  collect  the 
diffracted  light. 

In  addition  to  providing  a  proper  value  for  the  angular 
aperture,  all  particles  must  see  the  same  angular  aperture.  In 
other  words,  the  entrance  pupil  of  the  light  receiver  must  be 
large  enough  so  that  light  scattered  within  an  angle  equal  to  the 
acceptance  angle  of  the  receiver  will  be  detected  even  if  the 
light  come,  from  e  particle  at  the  outermost  edge  of  the  beam  and 
at  the  axial  location  farthest  from  the  receiver. 

A  system  employing  these  design  criteria  is  illustrated 
in  Figure  3. 

The  diameter  of  the  transmitter  tube  was  selected  as  a 
compromise  between  having  a  large  diameter  beam  to  sample  many 
particles  and  a  small  diameter  probe  to  minimize  flow  disturbances 
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when  inserting  the  probe  into  the  flow.  The  receiver  tube 
diameter  is  fixed  by  the  beam  diameter  and  the  required  entrance 
pupil  diameter.  The  angular  aperture  of  the  receiver  is  0.375 
degrees  while  the  theoretical  beam  spread  is  0,24  degrees. 

As  noted  above  the  problem  of  most  concern  to  the  proper 
interpretation  of  the  light  extinction  data  was  that  of  multiple 
scattering.  Prior  to  performing  the  full  scale  tests  in  the 
40"  x  18"  tunnel,  tests  were  run  in  a  12"  x  12"  tunnel  to  ascer¬ 
tain  if  multiple  scattering  was  important.  The  effects  of  mul¬ 
tiple  scattering  were  sought  by  comparing  average  densities 
measured  with  different  path  lengths.  The  measurements  required 
that  the  probes  be  inserted  into  the  flow.  These  small  scale 
tests  indicated  that  multiple  scattering  was  not  important  for 
particulate  densities  of  the  order  of  10/cm3  with  a  path  length 
of  26  cm.  The  full-scale  tests,  performed  at  higher  speeds  than 
the  small-scale  tests,  resulted  in  some  cases  in  a  much  larger 
cutoff  of  light.  Whether  or  not  the  Jouguer- Beer  law  holds 
under  all  full-scale  test  conditions  remains  to  be  established 
rigorously,  although  the  fact  that  the  full-scale  test  results 
are  consistent  among  themselves  (Section  3.3)  implies  that  the 
Bouguer-Beer  law  does  indeed  apply.  The  light  extinction  appa¬ 
ratus  was  built  to  allow  multiple  scattering  effects  to  be 
measured,  but  the  limitation  of  the  test  program  to  three  blows 
precluded  making  these  measurements.  Thus,  during  the  full 
scale  tests  the  length  of  the  light  beam  was  equal  to  the  full 
width  of  the  tunnel  -  18  inches.  Baseline  shift  during  low 
speed  runs  was  made  very  small  by  shielding  the  lenses  from 
the  dust- laden  particulate  flow,  but  at  higher  speeds  simple 
passive  shielding  was  les**  effective.  At  high  speeds  either  an 
air  flow  directed  outward  from  the  front  lenses  or  an  electro¬ 
static  screen  might  prove  helpful  in  reducing  baseline  shift. 


2.2.2  Wind  Velocity  Profile 

A  rake  containing  eight  pitot  and  three  static  tubes 
(Figures  4  and  5)  was  used  to  ireasure  the  velocity  profile  of  the 
wind.  The  pitot  tubes  employed  a  coaxial  design  (Figure  6)  to 
eliminate  clogging  of  the  tube  by  particulates.  A  fast  time 
response  was  designed  into  the  gauge  by  making  the  outer  passage¬ 
way  large  in  area.  By  placing  the  pressure  sensing  element  close 
to  the  front  of  the  gauge,  errors  in  the  measured  stagnation  pres¬ 
sure  arising  from  momentum  exchange  between  the  gas  and  particu¬ 
lates  have  been  eliminated. ^  Particulates  enter  into,  and  are 
decelerated  by,  the  inner  tube  without  affecting  the  measured 
value  of  the  stagnation  pressure.  A  bleed  hole  at  the  end  of 
the  inner  tube  allows  particulates  to  be  removed  from  the  inner 
tube;  thus,  the  gauge  test  time  is  extended  by  increasing  the  time 
to  fill  the  inner  tube  with  particulates.  The  bleed  hole  area  was 
made  equal  to  5%  of  the  probe  inlet  area  to  minimize  errors  in  the 
stagnation  pressure  measurement.6  A  pitot- static  tube  was  not 
built  although  a  triaxial  tube  can  be  readily  conceived.  The 
static  pressure  was  obtained  from  three  static  tubes  located  at 
different  positions  along  the  rake  which  also  held  the  eight  pitot 
tubes.  In  addition  to  supplying  a  static  reference  pressure  the 
three  static  tubes  were  used  to  determine  whether  or  not  stream¬ 
line  curvature  in  the  boundary  layer  Introduced  a  vertical  pressure 
gradient  into  the  flow. 

The  total  and  static  pressures  were  measured  with  Statham 
IW  131TC  differential  pressure  transducers  with  full  scale  ranges 
of  either  2.5  or  5.0  or  7,5  psla.  These  gauges  were  adequate 
for  the  required  task  but  required  a  careful  calibration  at  low 
pressures  employing  a  water  manometer.  The  total  differential 
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pressures  measured  during  the  tests  ranged  between  0.005  and  1.17 
psid.  The  acoustic  time  constant  of  the  pitot  tube  system  was  of 
the  order  of  7  msec,  which  was  much  faster  than  required  and  re¬ 
sulted  from  pre-test  underestimates  of  the  test  time. 


SECTION  3 


RESULTS 


3.1  VELOCITY  PROFILE  IN  THE  BOUNDARY  LAYER 

The  velocity  profiles  for  free  stream  velocities  of  111.8, 
237.8  and  376.3  fps  are  shown  in  Figures  9  to  11.  The  velocity 
within  the  boundary  layer  is  fitted  well  by  a  power  law  (Figure  12) , 
such  that  uay^a,  with  a  =  1.61,  1.36,  and  1.57  at  uq  =  111.8, 

237.8  and  376.3  fps  respectively.  A  semi- log  plot  of  the  velocity 
profile  (Figure  13)  reveals,  however,  that  the  flow  may  be  inter¬ 
preted  on  the  basis  of  the  law  for  flow  past  a  rough  wall.  The  three 

sets  of  data  intersect  at  y  «  1.65  in.,  u  =  30  fps.  This  intersec- 

9  10 

tion  is  qualitatively  similar  to  the  focus  of  Bagnold  and  Zingg 
but  occurs  at  a  greater  height  and  larger  velocity;  the  concept  of 
a  focus  is  therefore  approximate  at  best.  Moreover  the  constancy 
of  wind  speed  below  the  so-called  focus,  found  by  Bagnold  at  low 
free  stream  wind  speeds,  was  not  observed  at  the  high  free  stream 
wind  speeds  employed  in  these  tests  (Figure  13).  The  thickness  of 
the  velocity  boundary  layer  at  the  pitot  rake  station,  given  in 
Figure  13  by  the  intersection  of  the  two  straight  lines  that  charac¬ 
terize  each  velocity  profile,  varies  between  9.5  and  15.2  inches  as 
the  free  stream  velocity  varies  between  111.8  and  376,3  fps. 

The  shear  stress  velocity,  u*,  was  inferred  from  the  velocity 
profile  assuming  that  the  universal  constant  for  turbulent  flow  past 
a  wall,  5.75,  properly  described  the  experimental  data.  In  other 
words  u^  was  determined  by  applying  the  relation  du/d (log  y)  « 

5.75u*.  Thus  the  shear  stress  velocity  was  determined  to  be  16.3, 
41.0  and  61.0  fps  at  a  position  18  ft  from  the  leading  edge  of  the 
particulate  layer  at  free  stream  speeds  of  111.8,  237,8  and  376.3  fps 
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respectively.  That  the  constant  5.75  is  appropriate  is  also  shown 
by  plotting  u/ u^.  vs  log  y/(u^/2g)  (Figure  14).  Data  for  the  lowest 
speed  run,  111.8  f  ps ,  falls  on  the  straight  line  given  by  Owen^ 
as 

t  ■  5  75  logioir  +  9-7  <4> 

2g 

However  higher  speed  data  do  not  satisfy  Owen’s  relation.  Eq.  2 

12 

predicts  a  minimum  value  of  u/u^  =  5,  but  in  fact  much  of  the 
data  fall  below  this  value.  Boundary  layer  theory  (Appendix  I) 
makes  it  clear  that  Owen’s  relation  cannot  hold  near  the  leading 
edge  of  the  boundary  layer. 

3.2  RATE  OF  LOFTING  OF  PARTICULATES 

The  rate  at  which  particulates  are  lofted  by  a  high  speed 

.  _  C  /  A 

wind  is  given  experimentally  by  M  =  0.366*10  2  u  '  lbm/ (ft2-sec) 
(Figure  7).  The  variation  of  the  rate  of  mass  lofting  with  free 
stream  velocity  is  about  the  same  for  both  the  measured  and  com¬ 
puted  results,  which  implies  that  the  mass  lofting  rate  is  directly 
proportional  to  the  shear  stress  velocity  (Appendix  I).  The  rela- 
tion  between  M  and  u*  (Figure  8),  indicates  that  this  is  indeed  the 
case  with  a  least  squares  fit  yielding  M  -  (0.860*10-3  u^  -0.01) 
lbm/ (ft2-sec) .  The  implication  that  no  mass  is  lofted  for  shear 
stress  velocities  less  than  1.16  fps  is  in  qualitative  agreement 
with  threshold  values  determined  by  low  speed  tests. ^ 

That  the  measured  values  of  dM/di^are  the  same  as  those  com¬ 
puted  with  vertical  particle  velocity  equal  to  strongly  suggests 
that  the  initial  upward  velocity  of  the  particulates  is  proportional 
to  u^  throughout  the  range  of  wind  speeds  tested.  That  the  experi¬ 
mentally  determined  rate  of  lofting  is  an  order  of  magnitude  less 
than  the  computed  rate  implies  that  the  conversion  of  horizontal 
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impulse  at  the  ground  surface  to  vertical  momentum  of  particulates 
is  rather  inefficient.  It  should  be  borne  in  mind  that  the  param¬ 
eter  used  in  the  lofting  model  to  relate  the  initial  velocity  of 
the  particulate  to  the  shear  stress  velocity  was  extrapolated 
from  low  speed  tests  and  that  the  parameter  used  to  relate  shear 
stress  to  upward  momentum  was  chosen  to  maximize  the  lofting  rate. 
At  present,  it  appears  that  the  correct  value  of  the  latter  param- 

g 

eter  is  about  1/10  of  the  assumed  value. 

3.3  DENSITY  OF  LOFTED  PARTICULATES 

The  density  of  lofted  particulates  was  inferred  from 
light  extinction  measurements  (Figures  15  to  17),  using  single 
scattering  theory  (Section  2.2.1).  The  density  profiles  at  wind 
speeds  of  111.8,  237.8  and  371.3  fps  shown  in  Figure  18  correspond 
to  the  times  at  which  the  velocity  profiles  (Figures  9  to  11)  were 
measured.  A  linear  baseline  shift  correction  was  made  to  the  data 
in  the  calculation  of  the  densities,  but  the  rate  of  shift  of  the 
baseline  with  time  was  not  measured;  it  is  noted  that  the  error 
xn  the  derived  density  introduced  by  the  shift  is  made  less  severe 
by  the  logarithmic  relation  between  transmittance  and  particulate 
density.  The  error  does  not  exceed  25  percent. 

As  a  rough  check  on  the  accuracy  of  the  density  measure¬ 
ments  and  on  the  overall  consistency  of  all  of  the  measurements 
the  mass  flux  of  particulates  across  the  cross-section  of  the  tun¬ 
nel  was  calculated  and  compared  with  the  rate  of  removal  of  mass 
upstream  of  the  diagnostics  station  at  a  wind  speed  of  112  fps. 

The  flux  was  calculated  from  the  measured  densities  of  lofted 
particulates  (and  straight  line  extrapolations  to  y  -  0, 
n  ■  100/cra3  and  y  »  19.6",  n  *  0)  and  by  assuming  that  at  any 
given  height  the  particulates  were  transported  downstream  at  the 
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local  gas  velocity  (a  straight  line  extrapolation  of  the  log- 
log  velocity  plot  was  used  at  small  values  of  y) .  The  flux  cal¬ 
culated  in  this  manner  was  less  than,  but  within  15%  of  the 
lofting  rate  calculated  from  the  rate  of  depletion  of  the  partic 
ulate  bed. 
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SECTION  4 


CONCLUSIONS  AND  RECOMMENDATIONS 

* 

The  conclusion  of  greatest  practical  significance  is  that 
the  amount  of  wind- lofted  debris  is  only  1/10  the  maximum  value 
predicted  by  the  present  crude  theory.  That  aspect  of  the  model 
which  directly  relates  the  rate  of  lofting  of  the  mass  to  the  shear 
stress  velocity  is  borne  out  by  the  experiments.  The  gas  flow  in 
the  suspension  regime  is  qualitatively  similar  to  that  found  in  the 
saltation  regime  by  other  investigators;  namely,  the  velocity  pro¬ 
file  is  represented  well  by  the  law  for  flow  past  a  rough  wall. 
However,  Owen's  form  of  the  flow  lav;  was  found  experimentally  to 
hold  only  under  certain  conditions  In  addition,  the  design  phase 
of  this  program  made  it  appear  that  Owen's  relation  could  not  hold 
near  the  leading  edge  of  the  boundary  layer.  These  two  facts  point 
out  that  a  correct  calculation  of  particulate  lofting  must  account 
properly  for  the  growth  of  the  boundary  layer.  It  is  recommended 
that  any  future  study  of  lofting  include  an  analysis  of  the  two- 
phase  boundary  layer  growth  with  particulate  lofting.  Both  the 
skin  friction  coefficient  and  details  of  the  particulate  motion 
must  be  determined  to  characterize  adequately  the  properties  of  the 
boundary  layer.  The  photographic  determination  of  the  trajectories 
of  individual  particulates  is  made  difficult  by  the  large  concentra¬ 
tion  of  particulates  near  the  ground  surface  although  a  large  scale 
vortical  motion  of  the  particulates  is  revealed  by  high  speed  photo¬ 
graphs;  a  narrow  strip  of  particulates  may  make  photographic  analysis 
of  individual  particulates  possible.  If  so,  then  the  relationship 
between  flows  in  channels  of  different  widths  must  also  be  understood. 
A  high  frequency  particulate  counter  may  prove  useful  for  the  deter¬ 
mination  of  both  number  flux  and  particulate  velocity. 
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Transient  tests  similar  to  those  outlined  in  Appendix  I 
may  reveal  that  the  mass  of  material  lofted  behind  the  shock  wave 
is  greater  than  that  lofted  during  steady  flow, and  such  tests  should 
accordingly  be  undertaken.  Because  the  steady  flow  experiments  have 
now  been  performed,  at  least  in  part,  a  dual  purpose  tunnel  is  no 
longer  important  and  a  simple  shock  tube  would  suffice. 

Most  non-ideal  effects  such  as  a  compacted  surface,  moisture 
and  vegetation  would  make  the  lofting  rate  less  than  that  measured. 
Thermal  effects,  on  the  other  hand,  may  tend  to  increase  the  lofting 
rate.  The  effects  of  thermally-induced  shock  curvature  and  thermal 
gradients  should  be  investigated  to  complete  the  survey  of  lofting. 
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APPENDIX  I 


WIND  TUNNEL  DESIGN 
1.1  PRELIMINARY  DESIGN 

In  the  initial  design  of  the  lofting  facility  both  transient  and 
steady  state  tunnels  were  evaluated.  Important  factors  considered 
in  the  designs  were  the  rate  of  lofting  of  particulates,  the  inter¬ 
action  between  particulates  and  the  walls  of  the  tunnel,  the  thick¬ 
ness  of  the  boundary  layer  compared  to  the  height  of  the  tunnel,  the 
test  time  and  flow  stability.  The  transient  design  was  of  particular 
interest  at  the  outset  of  the  program  for  several  reasons.  For  example, 
by  proper  operation  of  a  blow-down  tunnel  the  flow  could  be  made  to 
pass  through  a  transient  phase  initiated  by  a  blast  wave,  and  then 
reach  a  period  of  quasi-steady  flow  limited  only  by  the  capacity  of 
the  blow-down  tanks  and  the  rate  of  depletion  of  the  particulate  bed. 
Thus,  the  blow-down  tunnel  would  not  only  supply  the  same  data  obtain¬ 
able  with  a  conventional  steady  state  tunnel,  but  would  also  furnish 
data  on  the  blast  wave  phase.  Moreover,  several  inviscid  calcu¬ 
lations  showed  that  the  blow- down  operation  would  reach  the  steady 
flow  in  the  shortest  possible  time;  namely,  steady  inviscid  flow 
would  be  achieved  directly  behind  the  blast  wave.  However,  the  in¬ 
clusion  of  viscous  effects  into  the  design  study  made  it  appear 
that  the  flow  would  be  unstable.  In  addition,  studies  of  the  two- 
phase  boundary  layer  made  it  appear  that  within  the  scope  of  extant 
theory,  a  high  speed  tunnel  had  to  be  very  large;  the  two-phase 
boundary  layer  was  estimated  to  be  much  thicker  than  an  air  boundary 
layer  at  the  same  free  stream  velocity.  Thus,  although  it  was 
thought  at  first  that  transient  tests  would  make  it  possible  to  in¬ 
vestigate  lofting  at  very  high  speeds  with  a  minimal  investment  in 
tunnel  facilities,  the  possibility  of  flow  instability  and  the  large 
tunnel  site  required  for  high  speed  tests  caused  the  transient  tun¬ 
nel  concept  to  be  discarded  in  favor  of  more  modest  speeds  in  a 
steady  state  tunnel. 
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1.1.1  Rate  of  Removal  of  Particulate  Mass 

The  scale  of  the  hardware  required  to  study  particulate 
lofting  at  a  given  wind  speed  is  governed  both  by  the  boundary 
layer  thickness  and  the  rate  of  removal  of  particulates  from  the 
test  bed.  Following  the  model  formulated  in  Reference  8  an  upper 
bound  to  the  particulate  mass  lofted  under  steady  state  conditions 
was  calculated  at  the  outset  of  the  program.  To  review  briefly 
the  lofting  model,  the  vertical  momentum  of  the  particulates 
lofted  from  a  loose  bed  was  set  equal  to  some  fraction  of  the 
horizontal  impulse  delivered  to  the  ground  plane  by  the  wind  in  a 
unit  time.  In  addition,  the  vertical  velocity  of  the  particles 
was  set  equal  to  some  fraction  of  a  quantity  u^,  known  as  the  shear 
stress  velocity.  For  the  purposes  of  estimating  wind  tunnel  loft¬ 
ing  both  fractions  were  set  equal  to  unity.  Under  the  two  stated 
assumptions  -  and  recalling  that  the  shear  stress  exerted  on 

the  ground  surface  by  a  wind  of  density  p  is  by  definition  equal 

8 

to  p  ul  -  we  find  that 
g  * 

Mvp  *  Mu*  -  Pgu3  At, 

where  v^  is  the  vertical  component  of  particle  velocity;  thus 

M  -  p  u*  At  (5) 

g  x 

A  proper  evaluation  of  Eq.  5  requires  that  the  relation  between 
the  shear  stress  velocity  and  the  free  stream  velocity  be  known. 

By  assuming  that  the  logarithmic  law  of  the  wall  holds  (Eq.  4), 
the  relation  between  uq  and  u*  can  be  obtained,  but  only  through 
its  dependence  on  the  boundary  layer  thickness.  In  fact  from 
Eq.  4  one  can  deduce  a  relation  between  u3/y  and  u/u*,  so  that 
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if  u  is  known  at  some  altitude  y  within  the  boundary  layer, 
then  u^  is  determined;  in  particular  u^  is  determined  explic¬ 
itly  given  tho  free  stream  velocity  uq  at  an  altitude  equal 
to  the  boundary  layer  thickness  6.  The  relation  between  u*/6 
and  uq/u^  (Figure  19,  with  y  =  6  and  u  =  uq)  indicates  that 
an  upper  bound  on  M  may  be  obtained  without  resorting  to  bound¬ 
ary  layer  calculations  by  setting  uQ/u^  =  5,  which  yields  the 
largest  value  of  u^  for  a  given  value  of  uq  (Section  3.2). 
Equation  (5)  then  yields  mass  lofting  rates  of  12.2  and  6.6 
lbm/ (fts-sec)  at  free  stream  velocities  of  800  and  430  fps, 
respectively.  If  particles  are  not  redeposited  on  the  test 
bed  then  the  rate  of  removal  of  mass  from  the  test  station  will 
be  equal  to  the  rate  of  lofting;  such  a  situation  is  encoun¬ 
tered  under  flow  conditions  that  characterize  the  suspension 
region  (Section  I)  which  were  in  fact  the  conditions  achieved  in 
the  experiments.  Thus,  maximum  rates  of  change  of  the  height  of 
the  test  bed  for  free  stream  wind  speeds  of  800  and  430  fps  are 
1.25  and  0.69  in. /sec,  respectively.  A  typical  test  bed,  which 
might  be  20  ft  long  x  1.5  ft  wide  would  then  experience  mass 
lofting  rates  of  0.275  and  0.148  tons/sec  at  800  and  430  fps, 
respectively.  Such  large  rates  of  mass  removal,  if  they  did 
indeed  occur,  would  render  Impractical  steady  state  tests 
at  high  wind  speeds  with  a  duration  of  a  second  or  more.  Skip¬ 
ping  ahead  for  a  moment,  the  results  of  a  more  detailed  calcu¬ 
lation  of  the  mass  lofting  rate  utilizing  boundary  layer  theory 
compare  closely  with  the  estimate  based  on  uq/u*,  but  as  shown 
by  the  experimental  results  both  mass  lofting  rates  are  gross 
overestimates. 
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1.1.2  Transient  Test 


The  large  mass  lofting  rates  obtained  as  discussed  in  Para. 
1.1.1  led  to  consideration  of  a  transient  test  with  a  limited  test 
time.  In  this  type  of  test  an  Initio],  blast  wave  is  followed 
by  a  uniform  flow.  The  use  of  a  blast  wave  is  advantageous 
not  onlv  because  aata  would  be  obtained  on  the  interaction  of 
such  a  wave  with  a  particulate  bed,  but  also  because  the  blast 
wave  will  lead  to  steady  flow  in  the  shortest  possible  time. 

Two  transient  flow  techniques  were  investigated  theo¬ 
retically  using  a  one- dimensional  variable- area  computer  code. 

The  first  scheme  utilised  a  test  section  with  an  initial  pres¬ 
sure  less  than  atmospheric.  The  blast  wave  strength  can  be 
adjusted  to  raise  the  test  section  pressure  to  1  atmosphere, 
thereby  matching  the  pressure  at  the  tunnel  exit  (the  acoustic 
impedances  would  be  slightly  different,  however,  depending  upon 
the  degree  of  shock  heating)  The  duration  of  the  test  is 
limited  either  by  depletion  of  the  sand  bed,  or  by  the  mass  flow 
characteristics  of  the  air  source.  The  results  of  a  typical 
machine  calculation  (Figure  20)  illustrate  the  idealized  sequence 
of  events  for  this  test.  The  plenum  chamber,  initially  at  a 
pressure  of  20.7  psia,  was  connected  to  the  test  section,  ini¬ 
tially  at  a  pressure  of  10,58  psia,  through  a  10°  nozzle  with 
A/ A*  -  6.25.  A  diaphragm  was  located  at  the  throat.  Upon  rup¬ 
turing  the  diapgragm  a  weak  shock  wave  moved  into  the  test  sec¬ 
tion  raising  its  pressure  to  13.52  psia  (the  particular  combina¬ 
tion  of  test  section  and  plenum  initial  pressures  did  not  result 


in  matched  exit  conditions).  At  the  same  time  a  rarefaction 
moved  back  into  the  reservoir,  inducing  the  gas  to  flow  into 
the  test  section.  After  1  millisecond  the  nozzle  became  over- 
expanded  causing  a  shock  to  stand  in  the  nozzle.  The  idealized 
inviscid  flow  at  a  given  station  in  the  test  section  reaches 
steady  flow  conditions  immediately  after  the  shock  passes  that 
station. 

The  second  scheme  utilized  a  test  section  at  an  initial 
pressure  of  1  atmosphere.  The  plenum  chamber  pressure  was  adjust¬ 
ed  to  achieve  the  desired  overpressure  and  flow  velocity  in  the 
test  section.  The  duration  of  the  test  was  limited  by  the  length 
of  duct  downstream  from  the  test  section  and  by  the  gas  sound 
speed.  Three  test  cases  were  calculated  to  investigate  the  ef¬ 
fects  of  nozzle  area  ratio  on  the  flow.  The  nozzle  area  ratios 
were  3.3,  4.3,  and  5.35.  The  pressures  within  the  flow  field, 
in  each  case,  are  plotted  in  Figure  21.  Each  case  was  designed 
to  yield  an  overpressure  of  approximately  10  psia  which  simu¬ 
lates  conditions  at  9500  feet  from  a  1  MT  blast.  As  in  the  pre¬ 
vious  case,  a  blast  wave  propagated  through  the  test  section  with 
a  stationary  shock  located  in  the  nozzle.  The  main  difference 
between  the  two  cases  lies  in  the  strength  of  the  blast  wave. 

The  results  of  six  calculations,  three  for  each  scheme,  are 
listed  in  Table  1.1. 
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TABLE  1.1 


SUMMARY  OF  TRANSIENT  TUNNEL  CALCULATIONS 


Velocity 

fps 

Static  Pressure 
psia 

Plenum  Pressure 
psia 

Initial  Pressure 
psia 

A/A* 

191 

13.5 

20.7 

10.58 

6.25 

413 

15.14 

46.7 

9.55 

6.25 

724 

11.75 

71.0 

4.78 

6.25 

431 

25.2 

82 

14.7 

5,35 

430 

25.2 

66 

14.7 

4.3 

416 

25.2 

50 

14.7 

3.3 

The  one- dimensional  model  does  not  give  a  complete  picture 
of  the  flow,  but  must  be  supplemented  by  a  knowledge  of  two- 
dimensional  effects.  In  particular,  the  ideal  discontinuity  stand¬ 
ing  in  the  nozzle  may  not  always  occur.  The  pressure  rise  pre¬ 
dicted  by  one- dimensional  theory  to  occur  across  a  discontinuity 
may  in  reality  occur  across  a  rather  large  distance  enclosing 

oblique  shocks  and  regions  of  mixed  subsonic  and  supersonic  flow. 

13 

This  extended  region  is  often  referred  to  as  a  pseudo- shock. 

The  pressure  rise  across  the  oblique  shock  system  often  encourages 

the  fluid  to  leave  the  walls  of  the  nozzle  with  resulting  flow 

separation.  Separation  in  steady  state  rocket  nozzles  has  been 

treated  in  the  literature^  The  maximum  theoretical  pressure 

ratio  that  can  be  taken  across  a  normal  shock  without  separation 
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is  approximately  1.67  for  a  flow  with  y  «  1.4.“  The  shock¬ 
boundary  layer  interaction  may  be  minimized  and  separation  pre¬ 
vented  by  removing  the  boundary  layer,  but  with  a  concomitant 


increase  in  experimental  complexity.  To  continue,  however,  even 

should  separation  occur  in  the  nozzle  it  is  almost  certain  that 

the  subsonic  flow  will  reattach  downstream  of  the  nozzle.  In 

fact,  certain  test  facilities  employ  a  choked  venturi  or  orifice 

plate  to  meter  flow,  relying  on  reattachment  downstream  to  smooth 
20 

out  the  flow. 

The  important  point  here  is  that  the  test  section  must  not 

be  too  close  to  the  nozzle.  Exactly  what  "too  close"  means  is  not 

quantitatively  known  at  the  present,  but  is  probably  of  the  order 

of  10  to  20  diameters.  In  addition,  a  limited  amount  of  experi- 
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mental  data  has  been  obtained  which  shows  a  normal  shock  stand¬ 
ing  for  approximately  50  millisec  in  a  highly  overexpanded  nozzle 
(Figure  22) .  These  data  were  taken  using  a  hot  shot  tunnel  employ¬ 
ing  a  hydrogen- oxygen- nitrogen  driver  and  a  nozzle  with  A/A*  =  1000. 
The  chamber  conditions  yielded  a  stagnation  pressure  of  approxi¬ 
mately  300  psia  and  a  stagnation  temperature  of  approximately 
6000°R.  The  downstream  section  was  evacuated  initially.  Although 
the  test  conditions  are  quite  different  from  those  needed  for  the 
lofting  test  these  data  suggest  that  a  standing  shock  may  occur 
for  a  limited  period  of  time. 

In  summary,  these  findings  indicate  that  during  a  transient 
test  a  normal  Lhock  may  stand  in  the  nozzle  under  some  flow  condi¬ 
tions.  In  a  steady  state  test  the  flow  will  almost  surely  separ¬ 
ate  if  the  pressure  rise  across  the  shock  is  greater  than  about 
1.67  ~  but  even  should  separation  occur,  the  flow  will  reattach 
downstream  of  the  nozzle. 

1.1.3  Boundary  Layer  Growth 

Initially  the  design  study  was  concerned  mainly  with  the 
problems  arising  from  the  large  rate  of  sand  lofting  anticipated 
in  a  high  spt.a  test.  The  problem  of  boundary  layer  growth  and 
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the  maintenance  of  an  undisturbed  free  stream  flow  with  minimal 
axial  and  transverse  pressure  gradients  was  also  undertaken.  The 
calculations  so  far  performed  apply  only  to  the  steady  state. 
Using  by  analogy  Walitt’s  results,  the  steady  state  calculation 
can  be  applied  to  a  transient  test  after  a  time  that  is  given 
approximately  by  2x/uq,  where  x  is  the  axial  distance  from  the 
leading  edge  of  the  sand  layer.  The  boundary  layer  thickness  is 
calculated  from 


dx  2  * 
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when  9  is  the  momentum  thickness  and  is  the  skin  friction  coef- 
icient,  with  the  assumption  that  the  velocity  profile  is  governed 
by  the  logarithmic  law  of  the  wall  (Eq.  4),  and  that  the  net  flux 
of  momentum  carried  by  the  sand  through  the  boundary  layer  is 
negligible 

Casting  Eq.  4  in  the  form 


(7) 


and  referring  to  Figure  19,  with  y  ■  6  and  u  *  u  ,  u  /u.  as  a  func- 

0  0* 

tion  of  (2g6/0)0/u*  may  be  approximated  by 


u. 


-0.87875 


log. 


r2g 


uo 


-  1.4985  log 


6  +  20.289 


(8) 


The  ratio  of  the  boundary  layer  thickness  to  the  momentum  thickness, 
6/0,  was  set  equal  to  (a+l)(a+2)/a  (Ref  23)  with  a  obtained  by  replot¬ 
ting  Zingg's10  data  as  a  log-log  relation.  The  average  value  of  a  was 


24 


23 

3.03.  Because  the  integral  equation  of  Von  Karman  was  employed 
to  determine  the  boundary  layer  thickness  it  is  felt  that  the  error 
introduced  by  setting  6/6  equal  to  a  constant  derived  from  a  power 
law  relation  is  small.  Introducing  En.  8  into  Eq.  6  and  integrat¬ 
ing  numerically,  the  relation  between  the  displacement  thickness, 

6*,  (Ref  23)  and  distance  along  the  particulate  bed  is  obtained 
for  various  values  of  the  free  stream  velocity  (Figure  23).  The 
calculation  predicts  a  finite  boundary  layer  thickness  at  the  lead¬ 
ing  edge  of  the  particulate  bed,  a  result  which  arises  from  use  of 
Eq.  4,  which  dictates  that  u/u*  ■  5  at  the  leading  edge  of  the 
boundary  layer.  Measurements  made  at  culmination  of  the  program 
show  u/u*  <5  at  some  locations  in  the  boundary  layer;  Eq.  4  can 
almost  surely  not  be  used  near  the  leading  edge  of  the  boundary 
layer. 

To  minimize  wall  effects  the  number  of  particulates  strik¬ 
ing  the  tunnel  roof  should  be  kept  as  small  as  possible.  The 
height  of  the  lofted  particulate  layer,  obtained  by  calculating 
particle  trajectories  in  the  boundary  layer,  is  plotted  in  Figure 
24  as  a  function  of  free  stream  velocity.  Figure  24  applies  to 
a  station  located  21  ft  from  the  leading  edge  of  a  particulate  bed 
containing  sand  grains  0.25  ran  in  diameter.  The  height  to  which  a 
particulate  is  lofted  increases  as  the  leading  edge  of  the  partic¬ 
ulate  bed  is  approached.  The  rate  of  decrease  of  lofted  height  with 
increasing  axial  position  decreases  as  the  wind  velocity  increases. 

At  a  low  velocity,  such  as  44  fps,  the  lofted  height  is  approximately 
h  the  value  of  u*/2g,  but  at  high  velocities,  such  as  350  to  430  fps, 
the  lofted  height  is  only  1/20  to  1/25  the  value  of  u*/2g. 

The  axial  increase  in  the  displacement  thickness  will  cause 
a  corresponding  increase  in  the  free  stream  velocity.  A  first  order 
approximation  for  the  increase  in  the  free  stream  velocity  relative 
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to  the  velocity  at  the  tunnel  entrance  for  various  tunnel  heights 
and  flow  velocities  is  plotted  in  Figures  25  and  26.  A  good  simu¬ 
lation  of  atmospheric  effects  requires  that  the  increase  in  the 
flow  velocity  at  any  axial  position  in  the  wind  tunnel  be  kept 
to  a  minimum.  However,  the  flow  velocity  increases  rapidly  with 
initial  velocity  for  a  fixed  tunnel  height;  whereas  for  uq  *  44  fps 
at  the  tunnel  inlet,  a  tunnel  18”  high  and  20'  long  will,  according 
to  the  calculation,  limit  the  increase  in  flow  velocity  to  about 
207;.,  a  tunnel  more  than  20  ft  high  is  needed  to  limit  the  increase 
in  flow  velocity  to  207o  when  uq  -  352  fps. 

The  great  disparity  between  the  relatively  low  lofting 
rates  determined  experimentally  and  the  upper  bound  provided  by 
the  theoretical  calculations  indicates  that  much  remains  to  be 
learned  al^out  boundary  layers  with  particulate  lofting. 
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APPENDIX  II 


TUNNEL  FLOW  CHARACTERISTICS 

The  relationship  between  the  time  for  the  tunnel  to  reach 
steady  state  gas  flow  conditions  and  the  rate  of  depletion  of  the 
particulate  bed  was  crucial  for  the  conduct  of  these  experiments. 
The  rise  time  of  the  tunnel  free  stream  velocity  was  fairly  long 
(Figure  27),  yet  the  relatively  low  lofting  rates  allowed  the  tests 
to  be  carried  out  without  modification  of  the  cross  wind  facility. 
The  test  time  for  the  purpose  of  calculating  lofting  rates  was 
determined  with  the  free  stream  velocity  arbitrarily  set  equal  to 
95%  of  maximum  free  stream  velocity,  a  convention  that  tended  to 
overestimate  the  lofting  rate. 
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Figure  1.  Lofting  Tunnel 
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MATERIAL:  304  STAINLESS  STEEL 
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Figure  6.  Coaxial  Pitot  Tube 
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Figure  7.  Lofted  Mass  at  Fixed  Position  in  Boundary  Layer  vs. 

Free  Stream  Wind  Speed.  Test  Bed  Composed  of  AFS  50-70 
Ottawa  Silica  Testing  Sand.  Mean  Diameter  of  Particulates: 
0.25  mm. 


37 


INCHtiS 


Figure  9.  Boundary  Layer  Particulate  Lofting  Teat  uq  -  111.8  fpa. 
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Figure  11.  Boundary  Layer  Particulate  Lofting  Test,  uq  .  376.3  fps. 
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Figure  20.  Calculated  Pressure  Distribution  in  Wind  Tunnel  at  Various  Times. 

C-as  Velocity  Behind  Shock  -  191  fps.  A/ A*  -  6.25.  Nozzle  Angle  * 
10  .  Initial  Pressure  in  Test  Section  »  1C. 58  psia. 
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Figure  21.  Calculated  Pressure  Distribution  In  Wind  Tunnel  at  10 

milllsec.  Nozzle.  Angle  «  10°.  Initial  Pressure  In  Test 
Section  -  14.7  pale. 
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(a)  Self  Luminous  Photograph  of  Flow  in  Convergent  -  Divergent 
Nozzle;  Luminosity  from  NaHCO^  Seed 


(b)  Karquardt  Nozzle  (Over-all  A/A*  =  1000) 


Figure  .’2.  Normal  Shock  in  Supersonic  Nozzle  of  Hot  Shot  Tunnel. 
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Figure  23.  Calculated  Displacement  Thickness 
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Figure  27.  Measured  Tunnel  Flow  Characteristics: 
Free  Stream  Velocity  vs.  Time. 


